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Abstract 

The form factors of the light pseudoscalar mesons are investigated in a dispersive 
formalism based on hadronic unitarity, analyticity and the OPE expansion of the 
QCD Green functions. We propose generalizations of the original mathematical 
techniques, suitable for including additional low energy information provided by 
experiment or Chiral Perturbation Theory (CHPT). The simultaneous treatment of 
the electroweak form factors of the 7r and K mesons allows us to test the consistency 
with QCD of a low energy CHPT theorem. By applying the formalism to the pion 
electromagnetic form factor, we derive quite strong constraints on the higher Taylor 
coefficients at zero momentum, using information about the phase and the modulus 
of the form factor along a part of the unitarity cut. 



1 Introduction 



Chiral perturbation theory |l] provides a systematic low energy expansion of the QCD 
Green functions in the sector of light mesons. Calculations up to two loops were done 
recently for the scattering amplitudes and the electroweak form factors in the chiral 
SU{2) x SU{2) limit §-§. A full 0(p e ) SU(3) x SU(3) calculation was also performed 
H for a specific combination of form factors which does not involve arbitrary renormal- 
ization constants. The result allowed to estimate the symmetry breaking corrections to a 
low energy theorem proposed by Sirlin |7|], which generalizes the Ademollo-Gatto theorem 

The fundamental properties of causality and unitarity are important ingredients in 
CHPT. Dispersion relations for scattering amplitudes and form factors or their inverses 
were used in performing low order calculations @], ||, or as a method of effective re- 
summation of the high order terms ||, [TIJ. In the present paper we shall consider an 



alternative dispersive method suitable especially for form factors [Tl|,[]12], which uses as 
input an information about the modulus along the unitarity cut, leading to constraints on 
the values inside the analyticity domain. The technique was for the first time combined 
with the OPE expansion of a QCD polarization function in | 13| . Recently, this formalism 



was successfully applied to the form factors describing the weak semileptonic decays of 
heavy mesons [fT4fl - pl|| . Quite strong constraints on the shape of these form factors near 
zero recoil have been obtained by combining the dispersive formalism with the predictions 
of the heavy quark effective theory (HQET) f22f , [p3fl . The predictive power of the tech- 
nique was considerably increased [XT], PU| , pif by the simultaneous treatment of several 
form factors with different unitarity thresholds, connected among them near zero recoil 
by HQET. 

It is of interest to apply this formalism to the form factors of the light pseudoscalar 
mesons, for which chiral symmetry predicts definite correlations for certain kinematical 
points. A first investigation of the weak form factors of the K — > tcIu decay in this 
framework was performed in [[U|. Similar techniques were applied in |24) in order to 



parametrize the modulus of the pion electromagnetic form factor in the time like region. 
In the present paper we perform an analysis of the electromagnetic and weak form factors 
of the 7r and K mesons, with emphasis on their low energy expansions. The purpose of the 
investigation is to see whether the dispersive approach is useful for testing the rigorous 
predictions of chiral symmetry and for constraining the free parameters of CHPT. To 
this end we present generalizations of the standard technique suitable for incorporating 
additional information provided by experiments or CHPT. 

The paper is organized as follows. In the next Section we review the dispersive for- 
malism: first we present the form factors of interest and some of their properties, then we 
describe the physical input of the method and the standard mathematical techniques used 
for optimally exploiting this input. In Section 3 we illustrate the simultaneous treatment 
of several form factors with different unitarity branch points by performing a test of a low 
energy theorem of CHPT. In Section 4 we derive constraints on the Taylor coefficients at 
zero momentum of the pion electromagnetic form factor. We first give the simple unitar- 
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ity bounds obtained from the standard formalism and then show how they are improved 
using additional information on the phase and the modulus of the form factor along a 
part of the unitarity cut. In the last Section we present some conclusions. 



2 Review of the dispersive formalism 
2.1 Definitions and notations 

We consider the electromagnetic form factors of the 7r and K mesons, defined as 

{K + {p')\Ji m \v + {p)) = {p + v')^{q 2 ) 
(K + (p')\j; lm \K+(p)) = (p + p'),F K+ (q 2 ) 

(K%p')\rJ m \K°(p)) = (p + p'),F K0 (q 2 ) , (1) 
and the form factors describing the weak semileptonic decay K — > ttIv: 
^\p')\J»\K+{p)) = (P + P%&\q 2 ) + (p-p')jl K \q 2 ) 

= L - U(q 2 ) + q»d nk (q 2 ) , (2) 



where q = p — p', P = p + p' and J^ m (J™) is the electromagnetic (weak) current: 

2 1 - 1 

J e J m = 3 ui^u - -dry ,4 - -s 7m s , J™ = s lyt u . (3) 

The form factors defined in (|1|) and (0) are analytic functions of real type in the complex 
plane t = q 2 cut along the real axis from the threshold of two particle production to 
infinity. The cut starts at t n = 4m^ for the electromagnetic form factors F n (t) , F K +(t) 
and F K o(t), and at t n K = {itl-k + ^k) 2 in the case of the weak form factors d W K(t) and 
fnK{t)- At t — the conservation of the electromagnetic current and the Ademollo-Gatto 
theorem || give 

F V (0) = 1, F K+ (0) = 1, F x o(0)=0, 
d nK (0) = (m 2 K - ml)f nK (0) , f nK (0) = 1 . (4) 

If we define the function 

A(t) = -i^t) + l -F K+ {t) + F K o(t) - f Kv (t) , (5) 
from (Hj) it follows that 

A(0) = 0. (6) 

This relation is valid up to terms quadratic in the chiral symmetry breaking parameters || . 
The Ademollo-Gatto theorem was extended by Sirlin M to values t ^ near the origin. 
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Sirlin's theorem requires in particular the vanishing of the derivatives of the function A(t) 
at t = 0: 

A'(0) = 0, A"(0) = 0,.... (7) 

In CHPT the renormalization constants cancel in the combination of the form factors 
entering Sirlin's function. A recent 0(p 6 ) calculation in full chiral SU(3) x SU(3) pertur- 
bation theory gave the value P 

A'(0) = -r| = (0.0033 ± 0.0005) fm 2 , (8) 

where r| is the "charge radius" of the Sirlin's form factor |J. This result shows explicitely 
the breaking symmetry correction to Sirlin theorem (0). 



2.2 Unitarity and dispersion inequalities 

Dispersive bounds on the above form factors are obtained by considering the vacuum 
polarization tensor : 



i / dxeW(0\T(Jt(x)J v {0))\0) = (q,q u - g^q 2 )U(q 2 ) + 9flu D(q 2 



(9) 



where J M denotes either the electromagnetic or the weak J™ current (the function 
D(q 2 ) vanishes in the electromagnetic case). From the asymptotic behaviour of QCD it 
follows that the derivative n'(g 2 ) of the amplitude IT(g 2 ) satisfies the dispersion relation 



n'(g* 



ImII(t + ie) 
(t-q 2 ) 2 



dt. 



(10) 



with the spectral function lmll(t+ie) given by hadronic unitarity. Using the definitions ([!]) 
of the electromagnetic form factors and taking into account the positivity lmll(t + ie) > 0, 
we obtain the inequality: 



lmll e i m (t + ie) > 



48tt 



t \ 3/2 

\F„{t)\ 2 9{t-t n ) 



48tt 



1 - — 

t 



\F K+ {t)\ 2 +\F K0 {t)\ 2 6{t-t K ) 



where t n = 4m 2 and tx = ^m 2 K are unitarity branch points. By inserting the inequality 
( p|) in the dispersion relation (|i~0D for an euclidian point q 2 = —Q 2 < we obtain: 



n; lm (-Q 2 )>^/ 



dt 



+ 



dt 



AM 2 J (t + Q 2 ) 2 



(t + Q 2 ) 2 
t K \ 3/2 



1 - 



4m; 



2\ 3/2 



t K 



1 - — 

t 



\F K+ (t)\ 2 + \F K0 (t)\' 



12) 
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A relation similar to ( |T2"[ ) can be written for the polarization function U W (Q 2 ) of the weak 
current. Using the definition (0) we obtain from unitarity and positivity 



ImII w (t + ie) > 



V 
48tt 

L ttK 



1 



I A 



+ \ 1/2 



IttK\ 



K 



t ( - ] 
L ttK 



9(t — Uk) , 



(13) 



where t^K = ( m K + "v) 2 , = (mj? — m^) 2 and 77 = | is an isospin factor fI5|| . By 
inserting (|13j) in a dispersion relation of the form fllCf ) we obtain 



00 



d/ 



X 



1 



T" 



1 - 



(t + Q 



2^2 



y- \ 1/2 



-n V 2 



t 



\fMt)? 



t 2 



(14) 



It is interesting to note that the inequalities ([12]) and ( fH|) can be added, leading to a 
single inequality 

.2 \ 3/2 



n'(-Q 2 ) > / 

K ^ J - 48tt 2 y 



487r 2 y (t + Q 



2\2 



1 - 



4m; 



+ 



+ 



dt 



487r 2 7 (t + Q 2 ) 2 



3/2 



dt 



48tt 2 y (t + Q 



2\2 



+ \ I/ 2 



i^ + (t)| 2 + |i^o(t)| 2 
( _ )X l/2 



t 



t 



X 



1 - 



t 



|/^(t)| a + 



t 2 



where 



(15) 
(16) 



n'(-Q 2 ) = n' eim (-Q 2 ) + K(-Q 2 ) ■ 

In the euclidian region Q 2 > 0, the functions U' elm (—Q 2 ) and n^(— Q 2 ) can be calculated 
by applying renormalization group improved perturbative QCD, with nonperturbative 
corrections included by means of operator product expansions (OPE). We used the ex- 
pressions given in 

n'(-Q 2 ) = 



1 i\ 1 



6 4/ vr 2 Q 2 
1 3\ 1 



71 



71 



- 77 + 



6 4/ 7r 2 g 4 



18 12/Q 6N tt 



1 , a, - 



+ [F A + 



10 



/3 2 7r 2 \ /a.(-Q 



2N\ 3' 



12 



7T 
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(^GG}+ - + 1 — m q (W), (17) 



q £ 



corresponding to rif = 3 flavours. The first and the second numerical constant in front 
of each term indicate the separate contribution of Tl' elm (— Q 2 ) and U' w (—Q 2 ), respec- 
tively. To evaluate these quantities we used the two-loop expressions of the running 
coupling a s (— Q 2 ) and the running mass m s {— Q 2 ) of the s-quark [P5[] , the perturba- 
tive parameters F 3 = 1.6398, F 4 = —10.2839, fli = —9/2 |25j and the condensates 
m q {^) = (0.200 GeV) 4 and (a s GG)/n = (0.45 GeV) 4 |2§. As concerns the value of 
the euclidian point Q 2 at which we calculate the polarization amplitude, it must be on 
one hand high enough to ensure the validity of the OPE expansion, and on the other hand 
small enough to provide a strong constraint on the form factors through the dispersion re- 
lation (|i0|) . In fact, the results turn out to be rather stable when Q 2 is varied in the range 
2 GeV 2 < Q 2 < 6 GeV 2 . The results reported in this work are obtained with the choice 
Q 2 = 2 GeV 2 , for which W elm (-Q 2 ) = 0.009546 GeV" 2 and W w (-Q 2 ) = 0.0133 GeV" 2 . 
With the l.h.s. of the relations (0), ( [14] ) or ([15]) known, these inequalities are integral 
conditions for the sum of the moduli squared of the corresponding form factors. 

It is useful to point out that an alternative integral condition for the electromagnetic 
form factors can be obtained by using as input, instead of OPE in the euclidian region, a 
lower bound on the hadronic muon anomaly. This method was first proposed in ||27| , and 
starts from the the vacuum polarization contribution of the hadronic part of the muon 
anomaly 



- / — lmU(t + ie)K(t) 
ix J t 



where 



K{t) = - [ - ^ "\ 9 dz . (19) 
w 7T J z 2 + (1 - z)t m 2 v ; 



z)t/ml 

By introducing (|T8|) in the unitarity relation (|ll|) we obtain the inequality 

,2\ 3 /2 




Ami 



\F n+ (t)\ 2 dt 



+lL J^(l-^) 3/ \\F K+ (t)\ 2 + \F K o(t)\ 2 ]dt, (20) 



t K 



which restricts the electromagnetic form factors along the unitarity cut. In the calculations 
we adopted the conservative lower bound aS£' > 7.5 x 10~ 8 [28], based on the experimental 



value, the weak radiative corrections and an estimate of the hadronic box diagram. 



2.3 Standard mathematical techniques 

In this subsection we briefly review the standard techniques for exploiting in the optimal 
way the conditions (|12[), (PI), 0I5[) or ( p0|) in order to obtain constraints on the form 
factors and their derivatives inside the analyticity domain. It is useful to introduce a 
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compact notation, denning the form factors Fi(t) , i = 1, .. , 5 as 



F x (t) = F n {t) , F 2 (t) = F K+ (t) , F 3 (t) = F K o(t) 
Fi{t) = U K (t), F 5 (t) = d« K (t) , 



(211 



and denoting by tj the unitarity branch points: 



ti — t n , t% — ix , ^3 



[22) 



We now apply to each integral with the lower limit ti appearing in the conditions ([12] 
(0)) (0) or ( P0|) the following conformal mapping 



(23) 



By this transformation the complex t plane is mapped onto the unit disk in the complex 
plane z, such that z(0) = and the unitarity cut t > t^ becomes the boundary \z\ = 1. 
For simplicity we denote Fi(z) = Fi(t(z)), where t(z) is the inverse of the mapping z(t) . 



Then the conditions (12) , fll4l), fllB]) or (|20[) can be written in the canonical form 



j 2tt 

— £ /K(C)F(C)| 2 d0<l, C = exp(z6 
27r i=i / 



(24) 



where / = 3 or / = 5. In Eq.([24]) W{(z) are analytic functions without zeros inside the unit 
disk, their modulus square on the boundary are equal to the weight functions appearing 
in front of the form factors, multiplied by the Jacobian of the conformal mapping (|23|) . In 
mathematical books these functions are called "outer functions" and are defined in terms 
of their modulus on the boundary as |29| 



wAz) 



exp 



, 2tt > 

1 r ( + z 



2tt J C-z 





]n\wi(C)\dO 



( = exp(U 



(25) 



In our case the functions Wi(z) have simple explicit expressions |13|-|pl||. For the global 
condition (|i5|) which involves all the form factors, the outer functions Wi(z) are 



'1 - <L 



WAZ 



16 V QrrtiW(-Q 2 ) (1 - zd { 



1,2,3, 



(26) 



w 4 (z) 
w 5 (z) 



32(1 - zJfl 2 V 6ttUII'{-Q 2 ) 



:i-rf.o 2 



[l + z) 2 VT=z~{l-zz-) 3 / 4 

(1 - zd A f 
- (l + z){l-zf 2 



32t 5 (l - z_fl 2 V 2irt 5 n>(-Q 2 ) (1 - zzJ) 1 /^! - zd 5 ) 2 ' 



(27) 
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where 



and 



y/U + Q' 2 -VU 
VU + Q 2 + VU 



tK 



1 ttK 



K 



For the alternative condition (M) the functions Wi(z) are replaced by : 



(2f 



(29) 



w i{ z) = - 



\ 



671^0^ 



(l + z) 2 VT^z~w K (z), 2 = 1,2,3, 



(30) 



where ty/c(;z) is an outer function defined as in (p5|) in terms of its modulus | w^:(C) I 
JC(C) on the boundary. 
The inequality (E¥) is a standard boundary condition in L 2 norm for the functions Fj( 



|29| , from which one obtains constraints on the size and shape of these functions at interior 
points. The results are expected to be stronger if one keeps in (^) the contributions of 
all the functions Fi(z) and exploits in addition the correlations among them provided by 
symmetries at some kinematical points. This technique proved to be very useful in the 
case of the weak form factors of heavy mesons, using HQET near zero recoil. In principle, 
a difficulty in applying the formalism to the excited states is the presence of the unphysical 
cuts below the unitarity branch points. In the case of heavy mesons, the unphysical cuts 
are well approximated by a few narrow resonances of known masses, which can be treated 
by adequate techniques with no assumption about the residua [M|-PT|. As we shall see 



in the next Section, the situation is more complicated in the case of the light mesons. 



3 Test of a CHPT low energy theorem 

In applying the relation ([24] ) we recall that the functions F\,F^ and F5, i.e. the elec- 
tromagnetic form factor of the pion and the weak form factors are analytic below their 
unitarity branch points, while the kaon electromagnetic form factors, denoted here as F 2 
and F 3 , have an unphysical cut along the region t n < t < fx, below the unitarity threshold 
of the KK production. In the variable z this cut is placed inside the unit disk \z\ < 1, 
along the segment — 1 < z < zq, where zq denotes the image of the two pion production 
threshold by the conformal transformation (|23| ) for the kaon form factors: 

Zq — — ^ — — — (31) 
y/tx — t w + \/Tk 

The dispersive formalism can be applied to this case only if some assumptions about the 
unphysical cut are adopted. For instance, if the phase of the kaon form factors along 
the unphysical cut is supposed to be known, the problem can be treated by means of 
a suitable modification of the standard techniques. Then one can use the inequalities 
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(|T2|) or ( [T5D to correlate the low energy expansions of the electromagnetic and weak form 
factors or to test some predictions of chiral symmetry. 

In this Section we will illustrate the procedure by testing the consistency of the rela- 
tions (^) proposed by Sirlin J7|, with the requirements of unitarity and perturbative QCD 
introduced in the condition (|T5|). We adopt for the phase of the kaon electromagnetic 
form factors a model similar to the one proposed in [Rtl for the pion form factor, and 



generalized to include a narrow to resonance, with the SU(3) relation C p = 3C U between 
the residua. We assumed therefore that the phases of the form factors F K + (t + ie) and 
F K o(t + ie) along the unphysical region are 



$f k+ (t) = arctg 



and 



$f k0 (t) = arctg 



2,m p T p (t) 
_A(m 2 p -t)\ ' "* 

3m p T p (t) 
2K " *) 



U < t < t 



K 



t n < t <C tfc 



(32) 



(33) 



We used the p width T p (t) given by the resonance chiral effective theory pl| , [p2 

,2 \ 3/2 



r,(t) 



m p t 
96^1 



1 



4m; 



(34) 



where /„- = 93.1MeV is the pion decay constant. 

We shall use now the information on the phase in order to remove the unphysical 
cut of the kaon form factors. To this end we shall use a so-called Omnes function p3| , 
which exactly compensates the known phase of the form factors along this part of the 
cut, leaving us with functions analytic below the unitarity threshold. In what follows 
we describe briefly the method. Since we treat all the form factors simultaneously it is 
convenient to introduce a collective notation, defining 



Oi(t) = exp 



Oi{t) = 
Si(t') 

t'(t' - 1) 



i = 1,4,5, 



dt' 



2,3. 



(35) 



where the functions 5i(t) are defined as S 2 (t) = Sp + (t) , 5 3 {t) = 6p^ (t) for t n < t < t K 



and are extended for t > t K as Lipschitz continuous functions |2{| . Of course, there is a 
large arbitrariness in such an extension, but the results are independent of the particular 
values of 5i(t) for t > tx- The reason is that, being analytic and without zeros in the 
t-plane cut along (t n , oo), and with regular values on the boundary, 02(f) and 0$(t) are 
outer functions, so by multiplying a class of functions with them, the class is not changed 
f29| . By construction, the phase of Oi(t+ie) coincides with 5i(t) . We also use the notation 
Oi(z) to denote these functions in terms of the variable z defined in (^). The functions 
Oi(z) are analytic in \z\ < 1, excepted for a cut along the segment (— 1,zq), where the 
phase of Oi(z — ie) is equal to 5i(t(z)) (we recall that the conformal transformation (|23| ) 
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maps the upper half t-plane onto the lower semidisk in the z-plane) . We define the set of 
functions fi(z) as 

f l {z) = F l {z)[O l {z)}~\ 2 = 1,.., 5, (36) 



where Fi(z) are the form factors defined in ([H]). From fl36|) it is easy to see that the com- 
plex phases of the form factors F 2 (z) and F$(z) along the cut (—1 , zq) are compensated 
by the phases of the Omnes functions Oi(z). Therefore, all the functions f%{z) are analytic 
of real type inside the unit disk \z\ < 1. By introducing fl3"6"|) in the L 2 norm condition 
(UM) we have 



2tt 



2tt 



i=l 



E / i^(c)o,(c)/,(c)i 2 ^ 



5 2n 

— £/ MCWC)/ 4 (OI 2 <w<i. C = ex P ( 

i=l n 



(37) 



where Wi(z) are the outer functions defined in (|26| ) and (27), and co>j(z) are additional 
outer functions which satisfy the relations 



k(C)l = |0i(C)|, C = exp(^6 



z = L...5. 



(3? 



They can be calculated in terms of their modulus on the boundary using the standard 
formula [29| given in (f25|), which is equivalent to 



Ui(t) = exp 



Vt^t r \n\Oi(t')\ 



-dt' 



, i = l,. ...5 



7T J y/p—fctf-t) 

It is convenient to introduce the new functions gi{z) as 

gi(z) = Wi(z)uJi(z)fi(z) = Vti(z)Fi(z) , i = 1, .. , 5 

where 

fti(z) = w^oj^z^O^z)]" 1 . 
Now, in terms of Qi{z) the inequality ( |5T| ) takes the canonical form 

-i 5 2n 

7T E/ki(C)| a <W<l, C = exp(^). 
27r i=i / 



(39) 



(40) 
(41) 

(42) 



By construction the functions <7i(z) are analytic inside the disk \z\ < 1 and can be ex- 
panded as power series 



n=0 



(43) 
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with real coefficients c n> i = c* i . By inserting these expansions in (021) we obtain 

5 oo 

EE<*<i- 

i=l 71=0 



(44) 



As we discussed in Section 2, the generalization of the Ademollo-Gatto theorem proposed 
by Sirlin requires the vanishing of the first derivatives at t = of the function A(t) defined 
in (|^), in the exact chiral limit. In order to test the consistency of this condition, let us 
put 

A {k) (0) = s k , k = l,....K, (45) 

where the derivatives are with respect to t, and Sk are given numbers of the order of 
magnitude of chiral symmetry breaking ||. Using ( f4"0"D we write A(£) as 



A(t) = J> [^{z)]' 1 gi {z) 



i=i 



where 



1 1 

2 2 

Therefore the derivatives entering (|45|) become 
A (fc) (°) = X>£CJ 

i=l n=0 

and finally the conditions fl4"5|) can be written as 

5 fc 



a 3 = 1 , a 4 = —1 , a 5 = . 



"(T^C*)]- 1 " 




"d fc -^(^)" 




dt n 


t=o 




t=0 



(46) 



(47) 



(4J 



Sk 



K. 



(49) 



j=l n=0 



where the numbers can be obtained in a straightforward way using Eq. (|48|), the 
definition (|4lD of the functions 0,%{z) and the connection (|23|) between the variables t and 
z. Of course, for each k only the Taylor coefficients c n ^ with n <k contribute to A( fc )(0). 

Our objective is to test the consistency of the relations (]45D with the inequality (|24]) 
provided by the dispersive formalism. The problem can be easily solved, since we ex- 
pressed both these conditions in terms of the Taylor coefficients c n> i, transforming in this 
way a functional problem into an algebraic one. We notice that the first coefficients 



Co,i = 9i{0) = ^(0)^(0) , % 



1,..,5 



(50) 



can be computed using the relations and pl|) . The remaining coefficients c n ^ n > 1 
are free. The relations ( |4"9"| ) can be viewed as a set of linear constraints for these 
coefficients, which satisfy in addition the quadratic condition (|4"4"D . A simple way to test 
the consistency of these relations for a given set of numbers s&, k — 1, .. , K is to define 
the quantity 



Mo 



mm 



i=l n=0 



(51) 
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where the minimization is with respect to the coefficients c^j, n > 1 , i = 1, ....5 which 
satisfy the set of linear conditions (f49|). It is clear that, if the minimal value /Xq is greater 
than one, the conditions ([44]) and ([49) can not be satisfied simultaneously, while if /ig is 
less than one the conditions are satisfied by at least one set of coefficients c n ^. A more 
detailed reasoning along these lines [34] shows that the inequality 

/io<l> ( 52 ) 

is a necessary and sufficient condition for the consistency of the relations ( f44"|) and fl49|). 

In order to solve the minimization problem (5"I) with the constraints fl4"9] ) we apply 
the technique of Lagrange multipliers. The Lagrangean of the problem is: 



5 5 oo K 

i=l i=l n=l k=l 



5 5 K 

£ ^fcO C 0,i + £ £ bkn C n,i - $k 



i=l 



i=l n=l 



(53) 



where k — 1, ..K are Lagrange multipliers. In (^) we indicated separately the coeffi- 
cients co,j which are known. Moreover, for convenience we extended the sum over n in the 
last term up to n = K, setting b^ n = for n > k. The Lagrangean £ is a convex function 
of the coefficients c n ^ its minimal value being given by the equations 



dC 

The solution of these equations is 



, n = 1, ...oo ; i 



L..5. 



(54) 



A" 



£A,6, 



kn i 



n < K 



1,..5 



k-l 







n > K + 1, i = 1, ..5 . 



(55) 



By introducing this solution in the conditions 
equations for the Lagrange multipliers 



we obtain the following system of 



K 



2_j Ukm^m — $k j k — 1, ...K . 



(56) 



m=l 



where the matrix U is defined as 



5 K 

TT —S^S^ h®h® 
u km / j / j u kj u mj 1 

i=l j=l 



k, m 



.K. 



(57) 



and the numbers Sk are related to the derivatives of the Sirlin function Sk defined in (|45"D 
by 



Sk = s k - £&lo c o,i , k = l,...K. 



(5f 



i=i 
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The system (|5B|) has the solution 



K 



h = E (U-^kmSm , k = l,...K, (59) 

m=l 

where U~ l is the inverse of the matrix U. By introducing these expressions in (|55|) we 
obtain the optimal coefficients c n> j for 1 < n < K as 

°n,i = E ^E( C/_1 )fe^ fo in) n = * = 1,-5. (60) 

m=l fe=l 

With these coefficients for n < K and c„ j; j = for n > K + 1 according to ( |55|) we 
obtain after a straightforward calculation the following compact form the solution /1q of 



the minimization problem (51): 



5 K 

^0 = E C li + E {U'^jkSjh , (61) 



i=l j,fc=l 

and the consistency condition ([52]) becomes 



E4 + ErV^<i- (62) 

i=l ^^=1 



We recall that in this inequality the coefficients co,j are known according to (|50D from the 
normalization of the form factors at t — 0, and the matrix {7 can be explicitely computed 
using fl5T|) and the numbers b^ n appearing in (^). We mention that these numbers include 
the physical input on unitarity and QCD , contained in the outer functions Wi(z), and the 
phase of the kaon form factors along the unphysical region entering the outer functions 
Oi(t) and u>i(t). As for the quantities Sk, they are related to the derivatives Sk of the 
Sirlin function, according to (0). With this input the inequality ((62|) is a constraint for 
the derivatives s& of the Sirlin function A(t) at t — 0, defined in d45f) . 

For illustration we evaluated the quantity //q for K = 3, using the value si = 0.0033 fm 2 
obtained recently in [§] and several choices for the higher derivatives S2 and S3. A few 
results are listed in Table 1. The sets of values for which /j,q is less than 1 are consistent 
with the dispersive bounds (in particular, this is true for S2 = S3 = and small deviations 
from these values, showing that Sirlin's theorem is verified within the present formalism). 
We give also some pairs (s 2 > S3) for which the quantity /j% is greater than one, and therefore 
the inequality (|62"D is violated. The results show an adequate consistency between the 
requirements of Sirlin's theorem on one hand, and QCD and analyticity on the other. 
However, in order to exploit the dispersive inequalities we had to adopt a model for the 
phase of the kaon form factors along the unphysical region. The presence of the unphysical 
cut reduces therefore the model independence of the conclusions. 
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s 2 (fm 4 ) 


s 3 (fm 6 ) 




0.0 


0.0 


0.04 


0.1 


-0.1 


0.69 


0.1 


0.1 


0.06 


-0.1 


0.1 


0.97 


-0.01 


0.3 


1.07 


0.05 


-0.3 


1.26 


-0.1 


0.5 


4.71 



Table 1: The quantity /Xq denned in fl6~Ij) for s± = 0.0033 fm 2 || and several values of 
the derivatives S2 and S3 of the Sirlin function. The values consistent with QCD and 
analyticity correspond to < 1. 

4 Bounds on the Taylor coefficients of the pion form 
factor at zero momentum 

In this Section we shall apply the dispersive formalism to the pion electromagnetic form 
factor F n (t). We consider the Taylor expansion of this function around the origin t = 

I 

F v (t) = 1 + lrlt + ct 2 + dt 3 + , (63) 

6 

where r\ is the radius of the charge distribution. In CHPT the calculation of the Taylor co- 
efficients requires the evaluation of higher pion loops, which introduce arbitrary renormal- 
ization constants §, §. One-loop CHPT § predicts c « 0.626 GeV" 4 , d « 2.30 GeV~ 6 . 
At two - loop level the coefficient c can not be calculated as it depends on an arbi- 
trary renormalization constant, and d = 4.1 GeV -6 ||. A fit of the ALEPH data 
p5| on the hadronic r decay rate with a Gounaris-Sakurai formula p6| (equivalent 
to the Pade version of the one-loop CHPT) gives c = 3.72 GeV -4 , d = 9.80 GeV" 6 . 
Similar values c = 3.9 GeV" 4 , d = 9.70 GeV~ 6 were obtained in [[37| by usual disper- 
sion relations. Other values proposed in the literature are: c = 4.1 GeV -4 [09(1 , and 



c = -7.5 GeV" 4 , d = 62.5 GeV" 6 @ 



In what follows we shall prove that the dispersive formalism imposes nontrivial con- 
straints on the Taylor expansion (|63"D, especially when combined with additional informa- 
tion about the form factor along a part of the time like region. For the present purpose 
we keep only the contribution of the pion form factor in the unitarity inequality (|TT|), ne- 
glecting the positive terms due to the K mesons. In this case the inequality ([24]) becomes 



1 2n 

^ / K(C)^(C)i 2 d0 < 1 , C = exp(i0) , (64) 
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where the outer function w w (z) is 



w (z) - {1 ~ d * f ! 1 (1 + zfVY ^~ Z (65) 



with d n = di defined in (P5[). If we use as an alternative input the lower bound on the 
hadronic contribution to the muon anomaly, the function w n (z) is 



w n (z) = — . 
16 



1 



(1 + zYVT^l w K (z) , (66) 



where Wk,{z) was defined below (|3~0|) . We consider now the function 

g(z) = w n (z)F 7T (z) , (67) 
which is real analytic in the unit disk \z\ < 1 and can be expanded as 

oo 

g{z) = Y J c n z\ (68) 

n=0 

with real coefficients c n = c * Then the condition (|64"D becomes 



E4<1- (69) 
n=0 



A constraint on the first iV coefficients appearing in the Taylor expansion (p3|) is imme- 
diately obtained by noticing that ( |69] ) implies 

N 

E4<1- (70) 

n=0 

The first N coefficients be expressed in a straightforward way in terms of the 

first N Taylor coefficients appearing in (|63|) . The explicit relations are obtained easily 
using the definition flST]) and the conformal mapping (p3|). For instance, c = ^^(0), 
ci = w' n (0) — 2rlt n w 7T (0)/3 etc, where the function w n (z) is defined in (|65|) and the 
derivatives are with respect to z. If we use as input the lower bound on the muon 
anomaly, the function w n (z) has the expression (|66|) . 

In particular, for N = 3 and a fixed value of the charge radius r\ the inequality ( ]7D| ) 
defines an allowed domain in the plane of the coefficients c and d appearing in the Taylor 
expansion (^). In Fig.l, the interior of the larger ellipse is the domain obtained using 
the standard value = 0.42 fm 2 and the QCD condition (|15|). This domain is very large, 
but one can see a certain correlation among the values of the Taylor coefficients. This 
feature becomes more stringent when higher derivatives of the form factor are taken into 
account. For comparison we indicated also in Fig. 1 the allowed domain obtained from 
the condition (p0|) on the muon anomaly. This domain is smaller, and in particular it 



excludes a pair of values for c and d proposed recently in . In the next subsections we 
shall improve the above bounds by implementing informations about the phase and the 
modulus of the pion form factor along a part of the unitarity cut. 
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Figure 1 : Allowed domains for the Taylor coefficients c and d of the pion electromagnetic 
form factor obtained from (fTo|). The large ellipse is obtained using as input the QCD 
expansion of the polarization function, the small one using the lower bound on the muon 
hadronic anomaly. 
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4.1 Improved bounds using the phase of the form factor along 
a part of the cut 



According to Watson theorem [EDI, along the elastic region t w < t < 16m^ the phase of 



the pion electromagnetic form factor coincides with the phase Sl(t) of the L — 1 , I — 1 
partial wave amplitude of the inr scattering. This amplitude was calculated in the frame 
of CHPT up to two loops |J. The expansion upon chiral loops is assumed to describe 
correctly the phase up to ~ 0.400 GeV. On the other hand, precise experimental data are 
available above ~ 0.600 GeV JO]]. New planned experiments on decay |4^[, ||43|| , [fPf 
will provide accurate information about the low energy tttt scattering amplitudes. 

It is of interest to incorporate in the dispersive formalism for the pion form factor 
the additional information on the phase. Unlike the situation encountered in Section 3, 
where the knowledge of the phase allowed us to remove the unphysical cut below the 
unitarity threshold, in the present case we must implement the phase along a part of the 
unitarity cut. We treat the problem adapting mathematical techniques used in ||45||-p8| . 
We assume that 

Arg[F^(t + ie)} = 5\{t) , U<t<t m , (71) 

where S\(t) is a known function and ti n denotes the threshold of the inelastic cut or a 
certain point up to which the phase is supposed to be known. As we shall discuss below, 
the technique can be easily adapted to more general situations, for instance when the 
phase is given along the range t\ <t < t in with t\ > t n , or along a region consisting of 
two disjoint parts, t n < t < t\ and t 2 < t < t in , with £ 2 > t±. 

We first express the condition ( |7T| ) in the variable z defined in (PB"|). To this end 
we recall that the branch point t n is mapped onto the point z — —1, and the upper 
edge of the unitarity cut in the t plane becomes the lower unit semicircle in the z plane. 
We denote by z in = exp(i9 in ) the image of the point t in — ie in the z-plane , and let 
T = {9 : 9 in < 9 < 2n — 9 in }. Then the condition ([H]) becomes 

limArg^re^)] = -81(9) , 9 in < 9 < vr , 

r— y1 

limArgfl^re^)] = 5\(9) , tt < 9 < 2vr - 9 in , (72) 

r — *1 

where we denoted 5{(9) = 5\(t(6)), with t(9) = t„ + t^cotg 2 (#/2) as follows from (H). 
We took into account the fact that the form factor is an analytic function of real type, 
which means that its phase is an odd function of 9, i.e. it satisfies the relation Sl(9) = 
— 8\(2-k — 9). It is useful to note that the same property holds also for the imaginary part 
of the form factor, while the real part and the modulus are even functions of 9. 

In what follows we shall derive the allowed domain for the Taylor coefficients of the 
expansion (|63"D, taking into account the inequality ( |64"D provided by the dispersive for- 
malism, and the additional relations ([721). As usual, the information about the phase is 
implemented by means of an Omnes function [53]. We can use the definition given in 
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(p5D, or equivalently we define 



O n (z) = exp 



2tt 



7T 



- / d9- 



1 — ze~ 



(73) 



where <SJ(0) is a Lipschitz continuous function such that 8\(9) = —8\{9) for 9 in < 9 < 
ii , 8\{9) — 8\{9) for n < 9 < n + 9 in , §1(9) being arbitrary outside T. We shall see below 
that the results are independent on the choice of S\(9) outside T. Using the Plemelj- 
Privalov relation |H9| 



lim — 

r-*l 7r 



2tt 



dO' 



no 1 ) 



1 _ re i(e-e>) 



F{9) 



2tt 



F{9') 



(74) 



where the last integral denotes the Principal Value, one can show that along the interval 
r the phase of the function coincides with the phase of the form factor. Therefore, 
by multiplying F n (z) with [0 7r (2)]~ 1 the phases compensate each other and the product 
is real along the elastic part of the cut. In the variable z this condition has the form 



Im lim 



i^(rO 



0, 9er. 



(75) 



By recalling the definition (|67l) 
the condition (|75|) in the form 



of the function g(z) and its power expansion (|68|) , we write 



Im lim 

r-+l 



wJre ie )OJre i6 



9{re w ) 



^c n Imlim [[lU(re i0 )]" 



1 r n inB 



n=0 



o, 9 e r 



where W(() is defined as 



w(0 = MC)o«(0> C 



(76) 



(77) 



The allowed domain of the Taylor coefficients of the pion form factor, which satisfy the 
conditions (^) and (|75|) can be found, like in Section 3, by means of an optimization 
problem. We first recall that the conditions (|64"D and (^) were written in the equivalent 
forms ( |69D and (|76| ) respectively, in terms of the Taylor coefficients c n . We consider then 
the quantity 



i-4 



oo 

min c 2 n 

n n=0 



(7E 



where the minimization is with respect to the coefficients c n which satisfy the condition 
(ff6|). As in Section 3 one can show that the inequality (0) is a necessary and sufficient 
condition for the consistency of the relations (^) and (|76|) . 
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We solve the constrained minimization problem ( [78] ) with the generalized Lagrange 
theory of optimization, based on Hahn-Banach theorem [50|. The Lagrangean is 



ra=0 



n=0 



£ = E 4 + - E c n Jim y A(0) |W(0) | Im [[W^)]" 1 rV" e J d9 , 



(79) 



where the function X(9) is a generalized Lagrange multiplyer |50| |. We can assume without 
losing generality that it is an odd function, since its product with another odd function 
(the imaginary part of an analytic function of real type) is integrated along a symmetric 
interval T. We denoted for simplicity W{9) = W(exp(i9)). The numerical factor in 
front of the integral in ([79|) and the modulus |W / (6 I )| inside the integral were introduced 
explicitely for convenience. 

We assume now that the first N coefficients c n have prescribed values, and perform the 
minimization of the Lagrangean with respect to the remaining coefficients c n , n > N + 1. 
The minimum condition 

n > N + 1 



dc n 







gives the optimal coefficients 

c n = --lim / X(9)\W(9)\lm \[W{9))- 1 r n e in " 

71 r— >1 J L 

r 

This expression can be written equivalently as 

c n = --lim / A(0) \W(6)\ [^(^l^^e 

7T r-»l J 



d9, n > N + 1 



0) 



(81) 



1 „n „—in9 



d9, n > N + 1 



(82) 



where we added a term which has a vanishing contribution by parity arguments. 

In order to find the Lagrange multiplier X(9) we introduce the optimal coefficients 
c n in the constraint (|76|). Recalling that the first N coefficients c n are assumed to have 
prescribed values, we write the condition ([75|) in the form 



N 

E c « Im 

n=0 









— Im lim 


_W(fl)_ 


r- >1 



|iy(y)|eW) 



W*( 



It is useful to write 
where from (|77|) it follows that 



7rW(0) 

1 

W{6) = |W(0)|e^, 

$(e) = ^) + 3J(fl). 



, ^T. (83) 

(84) 
(85) 

Here (f)(9) is the phase of the outer function w n (z), and 8{(9) is connected to the phase of 
the pion form factor as explained below (|73|) . By inserting (|^) in the relation (|83|) and 
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applying the Plemelj-Privalov relation fl75|) we obtain after a straightforward calculation 
the equation 

JV 

Y / c n sxn[n9-^(9))-X(9) 

n=0 
1 

+2^ 



\(ff) 



sin 






(9 - 9') - $(0) + $(0') 


sin 


'e-e r 

2 





de' = o , # g r 



(86) 



where the last integral is defined as the Principal Value. We obtained therefore a sin- 
gular integral equation for the Lagrange multiplier X(6). If the phase is Lipschitz 
continuous, the equation is of Fredholm type and can be solved by standard techniques. 
After solving this equation, the minimal norm /Xq can be computed by inserting in ([78]) 



the prescribed values c n for n < N and the optimal values (82) of c n for n > N + 1 
Taking into account the fact that c n = c* we obtain 



n 1 r r 

n\ = y cl + lim — lim / d9 / d0'A(0) X(9 

^ r-^X^r-^Xj J 



„ \W{9)\ \W{9')\ e KN+i)(e-e') 
W*(6) W*(6') 1 - re^-v') 



J7) 



By applying the Plemelj-Privalov relation (173) and using the integral equation fl86|) sat- 
isfied by the function A(0), we arrive finally at the expression 



N N 

»l = E 4 + E - / A (#) - *(*)]<w . 

n=0 n=0 71 



Using the fact that A(0) and $(0) are odd functions, the integral equation fl86| ) can be 
written in the form 



N 



Y^cn sin[n0 - $(0)] - A(0) 



n=0 



sin 


l( N + i) 


(0 - 9') - $(0) + $(0') 


sin 


6-0' 
2 





sin 


;(iv+i) (9+ 


9'- 


2tt) - 




- $(0')" 




sin 


'e+e 

2 


/ 

7T 








m 


= <\> 


(P)- 




9in < 



where 

Moreover, the expression ( jSli|) of //q can be written as 

AT AT 9 

n=0 n=0 



d0' = , 9 in <9<n, (89) 



(90) 



* N 2c r 
A*o = E 4 + E — / sinM - )]<W . 



(91) 
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We notice that the first term in Q9~H) represents the unconstrained minimum, obtained 
from ( [T0| ) if the coefficients c n for n > N + 1 were free. The second term in the expression 
of fiQ is positive and represents the improvement brought by the knowledge of the phase 
of the form factor along the region t n < t < t in . 

It is important to emphasize that only the values of along the interval V are 
required in the integral equation. Also, the modulus of the outer function does not 
appear (it was absorbed in the definition of the Lagrange multiplier A(#)), which means 
that the results are independent of the choice of the function 8\ outside the interval V . 

The equations (|89|) - ([H]) provide a simple numerical procedure for finding the allowed 
domain of the coefficients of the Taylor expansion (^): one starts with a set of values 
for the first iV Taylor coefficients. Using the relations (^)7]) and ( |68"D one finds the cor- 
responding coefficients c n , which enter the integral equation (j39|) . The solution \{6) of 
this equation is then used in fl9~l"|) to evaluate the quantity /Zq. Recalling that the allowed 
domain of the Taylor coefficients is described by the inequality fl52|) , the values taken as 
input are accepted or rejected if /i 2 , is less or greater than unity, respectively. 

From the above derivation the generalization to the case where the region V consists of 
two disjoint subintervals I\ and T2 is straightforward. This case is of interest when using 
the phase given by CHPT along T 1 = {t : U <t < (0.400 GeV) 2 } and the experimental 
data on T 2 = {t : (0.600 GeV) 2 < t < t in }. The resulting equations have the same form, 
the Lagrange multiplier being defined on Ti U T2- 

To illustrate how the bounds on the Taylor coefficients c and d are improved by the 
knowledge of the phase, we took the expression [BJJ 



6\{t) = arctg 



m p T p (t) 
m 2 — t 



(92) 



where T p (t) was defined in (0). At low energies above the threshold, this phase coincides 
with the one loop CHPT expression 

t ( Am 2 \ 3/2 

^-^rA 1 - J r) ' (93) 



while for t > (0.500GeV) 2 it is in very good agreement with the experimental data |41 
We assumed that the phase of the pion form factor coincides with (|92D along the region 
t n < t < t in , with ti n = 0.8 GeV 2 , which corresponds to 9i n = 0.6321. As input for the 
bounds we use the OPE expansion of the polarization amplitude, which led to the domain 
represented by the interior of the large ellipse in Fig.l. In this case the outer function 
w n (z) is given in ( |65|) and the function $(6*) defined in (|90|) has the expression 



$(0) = h arctg 



2di sinO -d\ sin 26 
. — d\ cos 6 + d\ cos 29 



-sl(e), e in <e<7r, (94) 



with d% defined in (P§|). In Fig. 2 the interior of the ellipse is the allowed domain in 
the plane (c,d), for r 2 = 0.42 fm 2 . Compared with the large ellipse in Fig.l one can see 
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20 40 

c (GeV -4 ) 

Figure 2: Improved domain for the Taylor coefficients c and d of the pion electromagnetic 
form factor, using the phase along a part of the boundary. 
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the considerable improvement brought by the knowledge of the phase along a part of 
the unitarity cut. A more realistic calculation using CHPT along an interval I\ and the 
experimental data on S\(t) along another interval is of interest and will provide precise 
model independent bounds on the Taylor coefficients. 

4.2 Improved bounds using information about the phase and 
the modulus along a part of the cut 

The modulus of the pion form factor along a part of the cut is known experimentally from 
the rate of e + e~ annihilation into pions |51| and the hadronic r decay []35|j . High precision 



data are available especially in the range 0.3 GeV 2 < t < 0.9 GeV 2 . In this subsection 
we derive the allowed domain of the Taylor coefficients by including some information 
about the modulus of the pion form factor along a part of the cut. 

It is a known mathematical fact that if the modulus and the phase of an analytic 
function are known exactly along a part of the boundary, then the function is in principle 
uniquely determined. Explicit formulas for recapturing an analytic function belonging 



to a certain class, from its restriction along a part of the boundary are available [52 



see also j|7| for explicit expressions and further references). These expressions are 



however very unstable numerically, reflecting the fact that the analytic continuation is an 
ill-posed problem []54| . More exactly, the formulas give the correct analytic continuation 



if the input values are known with infinite accuracy, but they lead to arbitrary predictions 
if these values are affected by errors. 

In what follows we shall show that even a nonoptimal use of the input information 
about the phase and the modulus leads to a considerable improvement of the bounds on 
the Taylor coefficients. Unlike the case considered in the previous subsection, where the 
phase could be given along several disjoint intervals, the method described below requires 
the knowledge of the phase along the whole range t n <t < t in . We start by writing the 
inequality (|12"D in the form 



1 7 dt / Ami \ 3/2 



48tt 2 7 (t + Q 2 ) 2 



Kit)] 2 < M{Q 2 ) , (95) 



where 



1 '?' dt / 4m^ 3/2 



Mm^-Q^-^jj^^-^} IF.WI'. (96) 

Assuming that the modulus is known along t n <t < t in , the quantity M.{Q 2 ) can 

be evaluated numerically. To be conservative, one can use a lower bound on the modulus, 
which does not spoil the inequality ([9*5]) and leads to a larger allowed domain for the 
Taylor coefficients. 

In order to incorporate the information about the phase we use the Omnes function 
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defined in ([73]), which can be written equivalently in the t variable as 



O n (t) = exp 



7T J t'it'-t) 



dt' 



(97) 



The properties of the function Sl(t) were discussed below ([73]), and we use it in order to 
remove the part of the cut below t — ti n . More precisely, writing as in ([36]) 



F*(t) = U{t)O n {t) 



(9* 



we see that the function f n (t) is real on the real axis below t in , since the phase of the form 
factor along t n < t < ti n is compensated by the phase of O n . Therefore f n (t) is analytic 
in the t-plane cut for t > ti n , and satisfies on the cut the condition 



dt 



Am 



2\ 3/2 



48ir 2 J (t + Q 2 ) 2 



\U(t)\ 2 \OAt)\ 2 < M(Q 2 ) , 



(99) 



which follows from fl9"5|) . But this condition can be simply exploited with the standard 
techniques described in subsection |2T3} We first map the t plane cut for t > t in onto the 
unit disk in the z-plane, using the conformal mapping 



z(t) 



\ftin t ytin 



(100) 



We define as in (P9|) the outer function 

co n (t) = exp 



y/t~=t f ln|O w (f)l 



7T 



-dt' 



V^tTn{t' - t) 

such that 1^(01 = 1^(^)1 for t > t in , and the outer function similar to 

q-d) 2 / r~ 



;ioi) 



with 



1 (l + z) 2 yT=I 

16 ^67iUM(Q 2 ) (1-zd^) 2 

\ftin + Q 2 - 



d 



y/tj/n. 

By introducing a new function g(z), defined as 

g(z) = w 1T (z)u 1T {z)f 1T (z) = w^u^p^z))" 1 F w (z) 
the inequality ( |9"9"D takes the canonical form 

1 2w 

— J \g(C)\ 2 d9 < 1 , C = exp(^). 



(102) 



(103) 



(104) 



(105) 
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Figure 3: Improved domain for the Taylor coefficients c and d of the pion electromagnetic 
form factor, using the phase and the modulus along a part of the boundary. 



By construction the function g(z) is analytic in the disk \z\ < 1 and can be expanded as 

oo 

~g(z) = Y / c n z n , (106) 

n=0 

where the coefficients c n satisfy the inequality 

oo 

E 3 < 1 , (107) 

n=0 

which follows from ( |105| ). The first N coefficients c n are connected in a straightfoward 
way through (|104j) to the first iV Taylor coefficients of the form factor F n (t) at t = 0. 
Therefore, the inequality 

N 

£ S n<l> (108) 
n=0 
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which follows from ( |107| ), defines an allowed domain in the plane of these Taylor co- 
efficients. In principle this domain can be further reduced: indeed, in addition to the 
boundary condition (|9^), the function f„ defined in fl98|) has known error-affected values 
along the region t n < t < ti n , where it is holomorphic. The inclusion of this additional 
information is not trivial, and we will not treat this problem here, using the simple in- 
equality (|108 ). In particular, for N = 3 and r\ fixed, this inequality gives the allowed 
domain of the parameters c and d of the expansion (|63"D. 

In a numerical application we used the same phase 5\ (t) as in the preceeding subsection 
and the modulus \F(t)\ from the ALEPH data |5|. Using Q 2 = 2 GeV 2 and t in = 
0.8 GeV 2 , we obtain for the quantity M{Q 2 ) defined in (||) the estimate M(Q 2 ) = 
0.00689GeV~ 2 . The result is not sensitive to the experimental uncertainties above the 
threshold t n , due to the phase space factor in the unitarity integral appearing in (|96|) . 
Reducing the value of \F(t)\ along t n < t < t in by 5% leads to an increase of M.{Q 2 ) by 
about 3% and to bounds on the coefficients c and d weaker by about 2%. The interior 
of the ellipse shown in Fig. 3 indicates the allowed domain obtained in this conservative 
situation. 

The combined information on the phase and the modulus is seen to restrict in an 
impressive way the values of the Taylor coefficients of the pion form factor. The allowed 
range of c is (0.25 , 7.57), and for each c the range of the values of d is very narrow. For 
instance, for c = 3.85 the parameter d is restricted to the interval (9.02, 11.30). The 
values c = —7.5 GeV -4 , d = 62.5 GeV -6 adopted in []38| are outside the allowed domain. 
We recall that these strong bounds are obtained with no specific parametrization and 
without any assumption about the high energy behaviour of the form factor. 



5 Conclusions 

In the present paper we investigated the form factors of the light pseudoscalar mesons in 
a dispersive formalism which uses as input the OPE expansion of the QCD polarization 
functions, combined with hadronic unitarity and analyticity. We derived constraints on 
the size and shape of the form factors, which are of interest in particular for testing the 
low energy predictions of CHPT. Generalizations of the original mathematical techniques, 
suitable for including additional informations about the form factors were developed. 

In Section 3 we performed a test of Sirlin's theorem 0, @, which requires the vanishing 
of a certain combination of form factors which is free of arbitrary renormalization constants 
in CHPT at t ^ 0. To this end we treated simultaneously in the dispersive formalism all 
the electroweak form factors of the 7r and K mesons. A difficulty is the large unphysical 
cut of the kaon electromagnetic form factors, which spoils the model independence of the 
results. This is in contrast with the case of heavy mesons, where the weak form factors of 
the ground state and of the excited ones have very close branch points and the unphysical 
cuts are well approximated by a few number of poles. In that case the simultaneous 
treatment of several form factors related among them by heavy quark symmetry led to 
strong model independent bounds near zero recoil fl6|| , [p0[| , pl| . In the present paper 
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we adopted for the phase of the kaon form factors along the unphysical cut a model 
inspired from the resonance chiral effective theory, which turned out to be consistent with 
QCD and Sirlin's theorem. This subsection has however mainly a pedagogical character, 
showing how to combine in an optimal way the dispersive formalism with the Omnes 
functions and the technique of Lagrange multipliers for the constraints at interior points. 

The main results are presented in Section 4, where we derived model independent 
constraints on the Taylor coefficients of the pion electromagnetic form factor. We showed 
how to implement in the dispersive formalism the knowledge of the phase of the form 
factor along a part of the unitarity cut, even when this part consists of several disjoint 
intervals. To this end we applied the Lagrange theory for functional optimization, which 
led to an integral equation for the generalized Lagrange multiplier. The knowledge of the 
phase considerably improves the simple bounds on the higher Taylor coefficients yielded 
by the dispersive formalism. In the present paper we illustrated this statement using 
a realistic model of the phase ||30|| , which reproduces the one-loop CHPT expression 



below 0.400 GeV and the present experimental data |4TJ above 0.500 GeV. The accurate 



data provided on the phase 5{ by the future K [4: experiments [|2|],j|3| will be a precious 
input to the formalism. In Section 4 we obtained also improved constraints on the Taylor 
coefficients if both the phase and the modulus of the form factor are given along a part 
of the cut, even if this information is used in an nonoptimal way. The interior of the 
ellipse shown in Fig. 3 represents the allowed domain of the coefficients c and d entering 
the Taylor expansion (|63|), obtained with the phase discussed above and the modulus 
from the hadronic r decay rate |35| along < t < 0.8 GeV 2 . Once this input informa- 



tion is adopted, the bounds are model independent, since they are not based on specific 
parametrizations and are free of any assumption about the high energy behaviour of the 
form factors. It is of interest to apply the techniques described in this paper to the higher 
Taylor coefficients. The results obtained so far suggest that strong correlations among 
these coefficients are expected, leading to restrictions on the arbitrary renormalization 
constants of CHPT. 
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